Elastic energy loss and longitudinal straggling of a hard jet 
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The elastic energy loss encountered by jets produced in deep-inelastic scattering (DIS) ofT a large 
nucleus is studied in the coUinear limit. In close analogy to the case of (non-radiative) transverse 
momentum broadening, which is dependent on the medium transport coefhcient q, a class of medium 
enhanced higher twist operators which contribute to the non-radiative loss of the forward light-cone 
momentum of the jet {q~ ) are identified and the leading correction in the limit of asymptotically high 
q~ is isolated. Based on these operator products, a new transport coefficient e is motivated which 
quantifies the energy loss per unit length encountered by the hard jet. These operator products 
are then computed, explicitly, in the case of a similar hard jet traversing a deconfined quark-gluon- 
plasma (QGP) in the hard-thermal-loop (HTL) approximation. This is followed by an evaluation 
of sub-leading contributions which are suppressed by the light-cone momentum q~ , which yields 
the longitudinal "straggling" i.e., a slight change in light cone momentum due to the Brownian 
propagation through a medium with a fiuctuating color field. 

PACS numbers: 12.38.Mh, ll.lO.Wx, 25.75.Dw 



> 
d 

oo 

o 



The modification of hard jets as they propagate 
through dense matter is now a rather sophisticated enter- 
prise [ll, H, H, 0, Hi J encompassing the study of jets prop- 
agating through cold nuclear matter in Deep-Inelastic 
scattering off large nuclei as well as jets in hot de- 
confined matter produced in the collision of heavy-ions 
at high energy Q. Unlike most approaches to jet mod- 
ification, the Higher- Twist expansion approach Q at- 
tempts to decompose the measured modification into a 
part computable using perturbative QCD (pQCD) and a 
set of non-perturbative "jet-transport-coefficients" which 
are then used to quantify the properties of the dense mat- 
ter. Since very few assumptions about the medium have 
been made, this formalism is equally applicable to both 
confined and deconfined matter. As a result, such a study 
has universal applications to a variety of heavy-ion ex- 
periments at the Relativistic Heavy-Ion Collider (RHIC) 
and the Large Hadron collider (LHC) and to DIS ex- 
periments at the Continuous Electron Beam Accelerator 
Facility (CEBAF) and the future Electron-Ion CoUider 
(EIC). To date, most jet modification calculations, in 
the higher twist scheme, have been limited to radiative 
energy loss and transverse broadening of the hard jet Q. 
Both these signals of jet modification are dependent on 
the transverse transport parameter q: which encodes the 
mean transverse momentum squared picked up by a hard 
jet per unit length of medium traversed 

In this Note, the longitudinal momentum loss of the 
jet due to elastic exchange with the medium will be ex- 
plored. We begin by first defining and clarifying the no- 
tion of elastic energy loss. It should be pointed out, that 
this is by far not the first attempt to explore this avenue 
of jet modification (see for instance Rcfs. (ol. ITol. [ill. [T3| ) . 
However, most of these computations assume a specific 
model of the medium: that of a thcrmalizcd plasma of 
quarks and gluons both within the HTL limit (as in 
Ref. [1, [13, [13]) and without invoking such limits [ll[. 



Our formulation, in this manuscript, assumes no partic- 
ular model of the medium and is thus related to that of 
Ref. fill . However, our results will differ from those of 
Ref. [llj l in that we identify a process which does not in- 
terfere with radiative energy loss and depends on a new 
transport coefficient. In so doing, the manuscript strives 
to attain a non-ambiguous definition of elastic energy loss 
within the various scenarios that may be encountered in 
current experiments on DIS and heavy-ion collisions. 

It is no surprise that the magnitude of elastic energy 
loss (and even radiative energy loss) is a frame depen- 
dent statement. While in heavy-ion collisions, the nat- 
ural frame of choice is the center of mass frame of the 
two colliding heavy-ions, there is more than one obvious 
choice in DIS on a large nucleus. As a result, in this 
Note, the focus will lie on the approximately boost in- 
variant quantity of fractional light-cone momentum loss, 
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where Eiost and pfg^f refer to the energy and z-component 
of the momentum lost by a hard parton with energy and 
z-component of the momentum given by E and p'^ . The 
qualifier, approximately boost invariant, is used as the 
above quantity is only invariant for the case of boosts 
restricted solely to the 3-direction, or the direction of the 
chosen component of momentum. Such a construction 
is somewhat alien to treatments of elastic energy loss in 
classical electrodynamics (See Chapter 13 of Ref. [l3l|). 
where the lost momentum is usually transverse to the 
direction of motion of the fast moving charge. 

Imagine that a hard quark is produced in the DIS of 
an electron on a large nucleus. Such an analysis is car- 
ried out with maximal ease in the Breit frame, where 
both the virtual photon and the l arg e nucleus approach 
each other at a large momentum [14l |. In the interest of 
simplicity, we imagine that the nucleus moves in the pos- 
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itive z-direction and has a large momentum where 
A is the atomic number of the nucleus. A hard quark 
with light cone momentum x bP^ is struck by the virtual 
photon which moves in the negative z-direction with mo- 
mentum 
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,0,0] , 
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where, the Bjorken variable has the usual value xb = 
Q"^ /{2p^q~) in terms of the virtuality of the virtual 
photon. In this frame, the struck quark has its four- 
momentum changed to I ~ XBP + q — [0, q~ , 0, 0], may be 
thought of as close to on-shell and moving in a direction 
opposite to that of the large nucleus. The differential 
hadronic tensor for the production of an on-shell hard 
quark with a momentum / in the interaction of a virtual 
photon on a large nucleus A is given as 



(PI I dl- 



= <VGx)<5(r-g-), 
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where, W^'^ represents the inclusive hadronic tensor at 
leading order and leading twist, i.e.. 



-Cp^g^Y^QlhixB) 
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where, we have used the light-cone notation —g^ 



g'^^Q^ ^ff'"'- The coefScient simply counts 



the number of nucleons in the nucleus that the jet may 
scatter off and is equal to A. 

In the remainder, we focus on the multiple soft reinter- 
actions of this quark with the soft glue field within the 
various nucleons. The weakly interacting picture of nu- 
cleons in the large nucleus (which is approximately valid 
at very high energies) is now imposed (see Ref. hj for 
details). Interactions of the quark with the various nu- 
cleons may henceforth be considered to be independent 
and thus uncorrelated. This step assumes a factoriza- 
tion of the production process of the hard quark and 
its later soft scattering off the glue field. As a result, a 
reader, more interested in the interactions of a quark with 
deconfined matter, may replace the individual nucleons 
with the prevalent degrees of freedom of the matter under 
study [il,tl6i]. 

Consider the soft rescattering of the hard quark with 
one such nucleon (or degree of freedom). In this limit of 
factorizcd interactions of the hard quark with the various 
nucleons, the effect of multiple interactions on the propa- 
gation of the hard quark may be iterated from the effect 
of a single interaction. This is not true in general and 
arguments will be forwarded which justify the use of this 
scheme in this particular case. By the effect on the hard 
quark, we specifically mean the effect on its three dimen- 
sional distribution in momentum space i.e., in terms of 
l~,l±. The aim is to identify the form of the distribu- 
tion (/)i(/~,/^), after the parton has traversed a certain 
length L given an initial distribution. 



The hadronic tensor with two gluon scatterings in the 
final state, after a few simplifications, may be expressed 
as 



M^'-=(-.9f)C,-,J dy,e 
X 9 



0. P+^O-(p|^(y-)l_^(0)|p) 



2 f dl dH^ _ + 2 dk d'^kjL 

(27r)35(r -q- - k-)S^{li_ - fc_i) Tr[ri''] 



2{q- + k- 



N, 



Tr 



7+7"Ur+r) + ..I^l_. -/^i 



X 7*^ I exp 



2(g- + 
k^ _ 

~'^2^(y ) + *y-L ■ ^ iy^k^ 



X (p2|^S(y-+y)4(r-)|p2). 



(6) 



In the limit of large energy, q~ — > oo, the distribution 
of k^ which is obtained by Fourier transforming the 
dependence of the product A°'{Y~ + y)A^(y~) will be 
dominated by values of k~ q~ . Under these limits, k~ 
may be dropped from the momentum dependent part of 
the integrand [third and fourth line of Eq. ([6])] and the 
integral over k~ may be performed using the phase factor 
exp (— iA;~j/+) which constrains the entire process to the 
negative lightcone. The resulting expression is then iden- 
tical to the well known expression obtained for the stan- 
dard treatment of transverse broadening with scattering 
off a single gluon in the final state The coefficient 
Cpp^, where \p) represents the nucleon with the struck 
quark and \p2) the nucleon which contains the soft glu- 
ons, accounts for A times the weak correlation between 
the two nucleons involved. This coefficient is derived in 
Ref. [3 and simplifications discussed in Ref. [T3 |. 

In the evaluation of transverse broadening in Ref. [TJ] , 
the two dimensional delta function S'^{1± — k±) was ex- 
panded in a Taylor series in k± and the coefficients of 
each term of the expansion became the gluon matrix ele- 
ments which appear in a product with the delta function. 
In the present case of non-radiative energy loss, we will 
follow a similar methodology and expand the delta func- 
tion S{1^ ~ q~ ~ k~) as a series in k~ and focus on the 
first derivative or the linear change in the distribution of 
l~ . Hence, the delta function may now be reexpressed as 



S{1- -q- -k-) = S{1- -q-) - 
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We use the equation above substitution for the 
(5-bmction in Eq. ([6|) and identify the coefficient (Ci) of 
the first derivative of the (5-function as the magnitude of 
the collisional energy loss. This may be schematically 
decomposed as. 



Ml',ll-)^S\ll_)5{l- -q-). 
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Ci =3C(fc)exp(fc.y)(y(2/)), 



(8) 
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i.e., a convolution of a purely momentum dependent part 
OC with the expectation of a purely position dependent 
piece y through a phase factor involving both. 

In the limit of <^ q~ and small coupling constant 
g, the momentum dependent piece 3C may be expanded 
in a series in k~ /q~ , i.e., 
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The piece which depends on the coordinates may be sim- 
ply expressed as V ~ A''°'{y)A^'^ {Q). In this expression, 
y = {y^ ,y~ ,y±) and assuming a weak dependence on 
, we have replaced Y~ with simply to save writing. 
Integrating over by parts, the overall factor of k~ in 
OC is converted into the derivative -~id~ acting on y. The 
terms in % are then ordered in powers of , the forward 
energy of the jet, which represents the largest scale in 
the problem. We begin by analyzing the leading term, 
which yields the coefficient of the first derivative of the 
longitudinal momentum (5-function as. 
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The 'oc' sign denotes that there is an over all factor that 
will arise from the decomposition of C^^^ in Eq. ([5]) [see 
Refs. [13, [Hi for details]. The coefficient (Ci) above 
may be further simplified by ignoring the small factor 
fc^/(2g~) in the exponent and integrating out the fc^ 
and the k~ to yield (5- functions over y^ and y\ . 

Unlike the case of transverse broadening [ij] , the two 
gluon matrix element is not manifestly gauge invariant. 
It may be cast in a gauge covariant form by noting that 
the gauge field at ^ oo is vanishingly small. In this 
limit, 



A'^+iY-) = / dz+d-A''+{Y-,z- 
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This yields a somewhat complicated expression for the 
elastic energy loss coefficient. In the limit of ^ oo, 
where one may ignore the k\/{2q~) in the exponent, the 
coefficient is given as, 
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The factor of /{2p'^), where, p^r is the nucleon density 
and p+ is the average forward momentum of the nucleon, 



is obtained from a factorization of the correlation coef- 
ficient Cpp^ ~ Cp X pjv/(2iP+), which is similar to the 
procedure used in Refs. [iJI to indicate the weak corre- 
lation between nucleons. 

In the derivation of the leading contribution to elastic 
energy loss in large nuclei, Eq. p2p is as far as one may 
proceed without invoking a model of the gluon distribu- 
tion within the nucleons. In the case of transverse broad- 
ening [l4| , a similar distribution involving the gauge field 
F~^-^ results; this is combined with the nucleon den- 
sity to define the transport coefficient q. In trans- 
verse broadening and energy loss in both large nuclei and 
heavy- ion collisions, q is often used as a parameter to fit 
with experimental data [l^ . We anticipate that a similar 
coefficient may also be motivated for the elastic energy 
loss, denoted as 



(13) 



2Nc 2p+ 

X J dy-dy+{p\ [tF--+{y-,0,0)F-'+{0,y+,0)] \p). 

The coefficient above represents the leading contribution 
to the elastic energy loss per unit length as encountered 
by a hard jet due to soft rescattering in the medium. 
There exist multiple other contributions such as energy 
loss due to hard rescattering, Compton scattering off 
hard gluons and subleading contributions (in q~ power 
counting) from Eq. ([9|. The first two contributions were 
calculated in Ref. [ll| and the emission of a hard for- 
ward parton in the final state leads to an interference of 
such contributions with radiative energy loss. As a result, 
these contributions are appropriately combined with and 
should be considered as a part of the radiative energy loss 
calculation [^. In this manuscript, the focus has been 
restricted only on those sub-leading contributions which 
do not have a hard on-shell gluon in the final state and 
thus do not interfere with radiative energy loss. Alterna- 
tively stated, we are only considering contributions where 
the parton loses longitudinal momentum q~ without pro- 
ducing another hard parton with large (— )-component of 
light cone momentum. 

As mentioned above, the effect of energy loss due to 
soft scattering off the glue field has been completely fac- 
torized from the production process of the hard jet. As a 
result, the coefficient e may be computed in any medium. 
In the case of elastic energy loss in a high temperature 
plasma, the coefficient in the equation above may be 
calculated exactly. Following the methods outlined in 
Ref. [l3|, we reabsorb the factor p/(2p+) into the defi- 
nition of a generalized medium state \n) and define the 
elastic loss coefficient as the expectation of the operator 
id^A^{y^)A~^{0) in the ensemble of thermal states as, 

ATraJdy-(n\e-'^"[id-A''+{y-)A''+{0)]\n) 
eHTL= . (14) 

Using the standard decomposition of two-point operators 
in finite temperature field theory [20| , we decompose the 
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correlator in the equation above, in the HTL hmit, as, 



ik+i 



2N, J (2^)4 



(15) 



where, is the Bose distribution function and p'^^^ rep- 
resents the (-|-+)-component of the gluon spectral density 
which is diagonal in color space, {pab = pSab)- 

The spectral density has the usual decomposition in 
terms of transverse and longitudinal components (using 
covariant gauge): 



PL 



(16) 



The projectors are further simplified on performing the 
y~ integral which yields (5(fc+). As a result, k^ = —2k^ 

and Pjt+ = = _p33 ^ 1^,2 1/1(^)2 _^ ;,2 1^ ugjj^g ^i^g 

standard notation (fc°/|fc| = x) [l^l, the expression for 
the elastic energy loss per unit length is given as. 



Ana, 
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X [{pT{\k\, k") - pL{\k\, fc")} {1 + riBik')} 



The equation above for the loss of light cone momen- 
tum l~ per unit of light cone path L~ is equivalent to 
the energy lost by near on-shell partons per unit length 
i.e. ., e ~ dE/dL. It may be evaluated numerically using 
the actual form of the gluon spectral density. In realistic 
calculations of the elastic energy loss at finite temper- 
ature, the form of the spectral density depends on the 
exchanged momentum of the gluon. For hard exchanges, 
one may simply evaluate the i-channel matrix elements 
for the scattering of a hard jet parton off a hard (^ T) 
medium parton. These are logarithmically divergent and 
are regulated at a soft scale p'^. Soft exchanges below 
this scale are appropriately calculated within the HTL 
formalism. On combining both these computations the 
arbitrary scale p^ is removed (see Refs. ^ [l^ for details 
and results of this procedure). 

The goal of this Note is not to rederive the results of 
Refs. [13 in their entirety but simply to show that a 
similar procedure may be carried out here with identical 
results. Adjusting for the choice of gauge and remov- 
ing the color factor Cp for the case of QED, Eq. (fT51) 
may be shown to be equivalent to Eq. (38) of the last 
article in Ref . Q . For the case of soft momentum trans- 
fers, one may use the known form of the HTL spectral 
densities [20| to compute e. Within this approximation, 
the upper limit of the |fc| integral is terminated at an 
appropriate scale of Qmax ^ T. With a choice of an 
as = 0.3 and a Debye mass of mjj ~ AnasT, this yields 
an e ~ 0.06 GeV/fm at a T = 300 MeV. In Refs. [a,[i3, 
the spectral density is rederived for the case of hard mo- 
mentum transfers; the difference between the results ob- 
tained from such a procedure and that from arbitrarily 



extending the HTL form of the spectral density to large 
momentum transfers is small (^ 20%) at large jet energy 
E (~ 40 - 50 GeV) [2l|. In this effort, we present this 
simplified estimate of the total elastic energy loss in a 
thermalized QCD medium by extending Qmax ~ V ET. 
This yields an energy dependent elastic energy loss per 
unit length as shown in Fig. [TJ The plot includes three 
different choices of T, and the coupling ag. The largest 
temperature and coupling yields the largest energy loss 
per unit length (e), with e dropping when the tempera- 
ture or the coupling is reduced. Previous efforts using the 
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FIG. 1: (Color Online) Elastic energy loss per unit length 
(dE/dL) as a function of the energy of the propagating parton 
for different fixed temperatures of the media and values of Qs . 
See text for details. 

higher twist formalism ignored this kind of mechanism of 
elastic energy loss. The current manuscript justifies the 
inclusion of such contributions within the basic formal- 
ism of higher twist elastic energy loss. Unlike the case of 
transverse broadening, there does not seem to be a simple 
classical analogue of this calculation, however, attempts 
in this direction are currently underway [l^. 

A closely related calculation, at the same order of ex- 
pansion in as, is the fluctuation of the elastic energy loss. 
This is obtained as the coefficient of the second derivative 
of the delta function of Eq. ([7]). Following a procedure 
almost identical to the steps carried out above, we obtain 
this coefficient in a thermal medium as, 



62 = 



Aira J dy 



-pii^pa — \- pa- 



2Nr 



(18) 



where, P° ^ has been defined in Eg. pT|) . Following a 
resummation procedure, as in Ref. |14l |. and by assuming 
that the correlation between F~^^ and P^ is vanishing, 
one may immediately postulate the diffusion equation for 
elastic energy loss as 



d(t){L-,l- 



dl- 



62- 



d{l-f 



, (19) 
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where the diffusion and drag coefRcients have been eval- 
uated above. 

In the interest of completeness, we will also evaluate 
the next-to-leading terms in Eq. ^ which provide addi- 
tional contributions to the clastic energy loss. These are 
however, further suppressed by the hard scale of q~ . The 
next-to-leading contribution to the elastic energy loss or 
expectation of k~ emanates from the contraction of the 
second term in 3C with ^ from Eq. ([9]). Evaluating this 
contribution, we obtain the coefficient of the first deriva- 
tive of the /" distribution as 

(X^y dY~dy~{p\d-A+iy-/2)d-A+{-y-/2)\p) 

^^JdY-dy-{p\F-+iy-/2)F-+i-y-/2)\p),i20) 

where, the standard shorthands of A" = and 
F = t'^F"' ^ have been used and a trace over color 
is implied. This term actually results in a slight gain in 
the longitudinal momentum fraction. In the derivation 
above, we have approximated that in the Breit frame 
(and in a covariant gauge) A'^ » A± ^ A~ and 
d~A~^ ~ F ''. While these contributions are suppressed 
by the large energy of the jet q~ , the operator products 
are not very different from those which constitute the 
leading contributions to elastic energy loss. As a result, 
such terms are only important at lower jet momenta. 



In this Note, we have presented an extension of the 
higher-twist expansion formalism of jet modification to 
include the effect of elastic energy loss. This was carried 
out by extending the formalism of transverse broadening 
in Rcf . [l^ , by generalizing the two dimensional distribu- 
tion of the propagating quark's transverse momentum to 
a three dimensional distribution including also its longi- 
tudinal momentum l~ . The equation governing the dis- 
tribution in includes both a diffusion term and a drag 
term. In contrast to other formalisms of jet modifica- 
tion, both these terms were evaluated at the operator 
level, independent of the details of the medium. The sole 
assumption used was that the color correlation length in 
the medium is small. The drag coefficient which yields 
the elastic energy loss per unit length was then evalu- 
ated in a thermal plasma in the HTL limit. The results 
obtained are consistent with similar calculations in other 
formalisms The final results, however, differed 

from previous attempts within the higher-twist formalism 
such as Ref. where elastic energy loss amplitudes in- 
terfered with radiative processes. In this Note, the focus 
has been solely restricted to those processes which do not 
interfere with radiative amplitudes. 

The author thanks U. Heinz, B. Miiller and G-Y. Qin 
for extensive discussions. This work was supported in 
part by the U. S. Department of Energy, under grant 
nos. DE-FG02-05ER41367 and DE-FG02-01ER41190. 



[1] M. Gyulassy and M. Plumer, Phys. Lett. B 243, 432 
(1990); X. N. Wang and M. Gyulassy, Phys. Rev. 
Lett. 68, 1480 (1992); R. Baier, Y. L. Dokshitzer, 
A. H. Mueller, S. Peigne and D. Schiff, Nucl. Phys. B 
483, 291 (1997); Nucl. Phys. B 484, 265 (1997); B. G. Za- 
kharov, JETP Lett. 65 (1997) 615; Phys. Atom. Nucl. 61 
(1998) 838 [Yad. Fiz. 61 (1998) 924]. 

[2] P. Arnold, G. D. Moore and L. G. Yaffe, JHEP 0206, 
030 (2002); S. Jeon and G. D. Moore, Phys. Rev. C 
71, 034901 (2005); S. Turbide, C. Gale, S. Jeon and 
G. D. Moore, Phys. Rev. C 72, 014906 (2005); 

[3] X. F. Guo and X. N. Wang, Phys. Rev. Lett. 85, 3591 
(2000); X. N. Wang and X. F. Guo, Nucl. Phys. A 696, 
788 (2001); B. W. Zhang and X.-N. Wang, Nucl. Phys. 
A 720, 429 (2003). 

[4] M. Gyulassy, P. Levai and L Vitev, Phys. Rev. Lett. 85, 
5535 (2000); Nucl. Phys. B 594, 371 (2001). 

[5] U. A. Wiedemann, Nucl. Phys. B 588, 303 (2000). 

[6] A. Airapetian et al. [HERMES Collaboration], Eur. 
Phys. J. C 20, 479 (2001); A. Airapetian et al. [HER- 
MES Collaboration], Phys. Lett. B 577, 37 (2003). 

[7] K. Adcox et al. [PHENIX Collaboration], Nucl. Phys. A 
757, 184 (2005); J. Adams et al. [STAR Collaboration], 
Nucl. Phys. A 757, 102 (2005). 

[8] A. Majumder, J. Phys. G 34, S377 (2007). 

[9] M. H. Thoma, Phys. Lett. B 273, 128 (1991). E. Braaten 
and M. H. Thoma, Phys. Rev. D 44, 2625 (1991); Phys. 
Rev. D 44, 1298 (1991). 
[10] M. Djordjevic, Phys. Rev. C 74, 064907 (2006). 



[11] X. N. Wang, Phys. Lett. B 650, 213 (2007). 

[12] G. Y. Qin, J. Ruppert, C. Gale, S. Jeon, G. D. Moore 
and M. G. Mustafa, Phys. Rev. Lett. 100, 072301 (2008). 

[13] J. D. Jackson, Classical Electrodynamics, John Wiley & 
Sons, Inc., New York, (1975). 

[14] A. Majumder and B. MuUer, Phys. Rev. C 77, 054903 
(2008); A. Majumder, R. J. Fries and B. MuUer, Phys. 
Rev. C 77, 065209 (2008). 

[15] E. V. Shuryak and I. Zahed, Phys. Rev. C 70, 021901 
(2004); Phys. Rev. D 70, 054507 (2004). 

[16] V. Koch, A. Majumder and J. Randrup, Phys. Rev. Lett. 
95, 182301 (2005); A. Majumder and B. MuUer, Phys. 
Rev. C 74, 054901 (2006). 

[17] X. f. Guo, Phys. Rev. D 58, 1 14033 (1998). 

[18] A. Majumder and X. N. Wang, arXiv:0806.2653 [nucl-th] . 

[19] A. Majumder, C. Nonaka and S. A. Bass, Phys. Rev. 
C 76, 041902 (2007); S. A. Bass, C. Gale, A. Ma- 
jumder, C. Nonaka, G. Y. Qin, T. Renk and J. Ruppert, 
larXiv:0808.0908l [nucl-th] . 

[20] J. I. Kapusta and C. Gale, Finite-temperature field the- 
ory: Principles and applications, Cambridge, UK: Univ. 
Pr. (2006) . 

[21] G.-Y. Qin, Ph.D. Thesis. 

[22] R. B. Neufel d, B. MuUer and J. Ruppert, 
larXiv:0802.2254l [hep-ph]; R. B. Neufeld, 
larXiv:0805.0"385] [hep-ph]. 

[23] M. Luo, J. w. Qiu and G. Sterman, Phys. Lett. B 279, 
377 (1992); Phys. Rev. D 50, 1951 (1994). 



